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Abstract — Head losses are very common in penstock trifurcations. In this paper, six cases of 3D asymmetric trifurcations 
have been modeled with main pipe length & diameter of 1.3716m & 0.0254m, respectively and branch pipe lengths & 
diameters of 0.762m & 0.0196m, respectively. Volumetric flow rates, velocity magnitudes, dynamic and total pressure 
contours and their values have been computed. Energy loss coefficients have been computed for branch pipes for an input air 
velocity of 3m/ s by pressure data obtained from the CFD analysis. The maximum values of velocity magnitude, dynamic and 
total pressures are observed in the branch-2 and head losses in branch-2 are relatively less. 

Keywords — Head losses , Energy loss coefficients, 3D asymmetric trifurcations. 

I. Introduction 

In penstocks used for hydropower projects, trifurcations along with the other components, help in producing electricity. 
These trifurcations supplement water supply to multiple turbines at the same time. Despite having the economical advantage 
over independent systems, even this system is not free from losses. The comparison of velocity magnitudes, dynamic and 
total pressure contours and determination of head loss coefficients in the branch pipes sums up the interest of this study. 

II. Domain 

A total of six cases have been modeled using Gambit 2.4.6, with each model displaying an asymmetry about the central 
branch axis. The dimensions used for the six cases in the current analysis are shown in the table 1. 


Table 1 

Dimensions of Asymmetric Trifurcations 


Dimensions 

In Meter 

Length of Main Pipe (L) 

1.3716 

Diameter of Main Pipe (D) 

0.0254 

Lengths of branch pipes (li, 1 2 , 13) 

0.762 

Diameters of branch pipes (di, d 2 , d 3 ) 

0.0196 


The angle between the branch-2 and branch- 1 is termed “a” and the angle between the branch-2 and branch-3 is termed “p”. 
The angles used for the six cases are shown in the table 2. 


Table 2 

Angles “a” and “P” for the Six Cases 


Case No. 

1 

2 

3 

4 

5 

6 

Angle “a” in Degrees 

5 

10 

15 

10 

15 

20 

Angle “P” in Degrees 

10 

15 

20 

5 

10 

15 


A model created using Gambit for the case 2 of this study is shown in figure 1 . 

Mesh of T-Grid type using Tet/Hybrid elements was generated for all the six cases in Gambit. A typical mesh generated for 
the case 3 of this study is shown in figure 2 
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Fig. 1: Gambit Model of Trifurcation Fig. 2: Meshed Trifurcation 


The mesh sizes for all the cases are enlisted in table 3. 


Table 3 

Mesh Sizes for the Six Cases 


Case 

Mesh Size 

No. 

Cells 

Faces 

Nodes 

1 

469301 

978796 

99471 

2 

320336 

674755 

71535 

3 

320336 

674755 

71535 

4 

316436 

666226 

70469 

5 

321042 

676181 

71632 

6 

314140 

662439 

70526 


The boundary conditions for all the six cases were applied as shown in the figure 3 and table 4. 





Fig. 3: Application of the Boundary Conditions at Inlets and Outlets 

Table 4 

Boundary Conditions For The Domain 


Pipe 

Entity 

Boundary Condition Type 

Magnitude 

Main Pipe 

Face 

Velocity Inlet 

3 m/s 

Branch Pipe-1 

Face 

Pressure Outlet 

OPa 

Branch Pipe-2 

Face 

Pressure Outlet 

OPa 

Branch Pipe-3 

Face 

Pressure Outlet 

OPa 


Ansys fluent 13.0 was the solver used for the analysis. The details of the fluid properties and solver parameters are given in 
the tables 5 and 6 respectively. 
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Table 5 


Fluid Properties 


Fluid 

Fluid Properties 

Density 

Viscosity 

Air 

1.225kg/m 3 

1.7874 x 10' 5 kg/m-s 


Table 6 


Solver Parameters, Initialisation and Calculation Details 


Pressure- Velocity Coupling 

Simple 

Spatial Discretization 

Gradient 

Least Square Cells Based 

Pressure 

Second Order 

Momentum 

Second Order Upwind 

Solution Initialization 

Initialisation Methods 

Standard Initialisation 

Reference Frames 

Relative to Cell Zone 

No. of Iterations 

1000 


The analysis has been carried out for the six cases under the following assumptions: 

• No slip condition; which means that the relative velocity of the fluid at the solid boundaries is zero. 

• The fluid flow is incompressible. 

• Air is a Newtonian fluid. 

• Steady flow occurs. 

III. Volumetric Flow Rates, Velocities and Dynamic & Total Pressure Values 

The values of velocity magnitudes, dynamic and total pressures have been obtained for the branched flow as well as inlet for 
the surface areas of all the six cases from fluent. These values are tabulated in the tables 7-12. 


Table 7 

Volumetric Flow Rates, Velocity Magnitudes, Dynamic and Total Pressures for Case 1 


Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1.5 x 10' 3 

3 

5.3621 

15.7877 

Branch Pipe-1 

4.32 x 1 0 4 

1.44 

1.4374 

1.4466 

Branch Pipe-2 

6.94 x 10 4 

2.32 

3.6253 

3.6428 

Branch Pipe-3 

3.82 x 10 -4 

1.28 

1.1392 

1.1464 

Volumetric Flow Rates, Velocity M 

Table 8 

[AGNITUDES, DYNAMIC AND TOTAL PRESSURES FOR CASE 2 

Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1 .5 x 1 0 

3 

5.3892 

16.0206 

Branch Pipe-1 

4.4 x 10' 4 

1.47 

1.4711 

1.4851 

Branch Pipe-2 

7.05 x 10' 4 ' 

2.36 

3.6804 

3.7145 

Branch Pipe-3 

3.61 x 10 -4 

1.21 

1.0035 

1.0092 

Volumetric Flow Rates, Velocity M 

Table 9 

[AGNITUDES, DYNAMIC AND TOTAL PRESSURES FOR CASE 3 

Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1.5 x 10' 3 

3 

5.3957 

16.2117 

Branch Pipe-1 

4.2 x 10' 4 

1.41 

1.3503 

1.3603 

Branch Pipe-2 

7.02 x 10' 4 

2.35 

3.6775 

3.7062 

Branch Pipe-3 

3.84 x 10 -4 

1.29 

1.1277 

1.1363 
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Table 10 

Volumetric Flow Rates, Velocity Magnitudes, Dynamic and Total Pressures for Case 4 


Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1.5 x 10 -3 

3 

5.3956 

15.7296 

Branch Pipe-1 

3.6 x 10' 4 

1.21 

1.0018 

1.0100 

Branch Pipe-2 

6.66 x 10~ 4 

2.23 

3.3027 

3.3298 

Branch Pipe-3 

4.8 x 10' 4 

1.61 

1.7367 

1.7520 

Table 11 

Volumetric Flow Rates, Velocity Magnitudes, Dynamic and Total Pressures for Case 5 

Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1.5 x 10 -3 

3 

5.3892 

16.0252 

Branch Pipe-1 

3.84 x 10' 4 

1.29 

1.1389 

1.1488 

Branch Pipe-2 

6.89 x 1(T 4 

2.31 

3.5319 

3.5633 

Branch Pipe-3 

4.32 x 10' 4 

1.45 

1.4160 

1.4296 

Volumetric Flow Rates, Velocity M 

Table 12 

[AGNITUDES, DYNAMIC AND TOTAL PRESSURES FOR CASE 6 

Surface Area 

Volumetric Flow Rate 

(m 3 /s) 

Velocity Magnitude 

(m/s) 

Dynamic Pressure 
(Pa) 

Total Pressure 
(Pa) 

Inlet 

1.5 x 10 -3 

3 

5.3890 

16.2039 

Branch Pipe-1 

3.95 x 1(T 4 

1.33 

1.1957 

1.2054 

Branch Pipe-2 

6.86 x 10' 4 

2.3 

3.4997 

3.5283 

Branch Pipe-3 

4.25 x 10 -4 

1.42 

1.3707 

1.3799 


IV. Velocity Magnitude Contours 


■ 4 56e+00 
4 34e+00 
4 11e+00 
3 88e+00 


SYS 


3 65e+00 
3 42e+00 
3 19e+00 


2 97e+00 



Contours of Velocity Magnitude (m/s) 

May 16. 2016 


ANSYS FLUENT 13.0 (3d. pbns. ske) 



Contours of Velocity Magnitude (m/S) 

May 16.2016 


ANSYS FLUENT 13.0 (3d. pbns. lam) 


Fig. 4 Velocity Magnitude Contour-Case 1 (m/s) 


Fig. 5: Velocity Magnitude Contour-Case 2 (m/s) 
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Contours of Velocity Magnitude (m/s) 

May 16.2016 


ANSYS FLUENT 13.0 (3d. pbns. lam) 


Contours of Velocity Magnitude (m/s) 

May 16.2016 


ANSYS FLUENT 13.0 (3d, pbns. lam) 


Fig. 6: Velocity Magnitude Contour-Case 3 (m/s) 


Fig. 7: Velocity Magnitude Contour-Case 4 (m/s) 



Fig. 8: Velocity Magnitude Contour-Case 5 (m/s) 


Fig. 9: Velocity Magnitude Contour-Case 6 (m/s) 


V. Dynamic Pressure Contours 



Contours of Dynamic Pressure (pascal) 


May 15.2016 
ANSYS FLUENT 13.0 (3d. pbns. lam) 



Contours of Dynamic Pressure (pascal) 


May 16. 2016 
ANSYS FLUENT 13.0 (3d. pbns. lam) 


Fig. 10: Dynamic Pressure Contour-Case 1 (Pa) 


Fig. 11: Dynamic Pressure Contour-Case 2 (Pa) 
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Contours of Dynamic Pressure (pascal) 


Jun 25,2016 
ANSYS FLUENT 13.0 (3d. pbns, lam) 
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Contours of Dynamic Pressure (pascal) 

May 16,2016 


ANSYS FLUENT 13.0 (3d. pbns. lam) 


Fig. 13: Dynamic Pressure Contour-Case 4 (Pa) 



Contours of Dynamic Pressure (pascal) Jun 25. 201 6 

ANSYS FLUENT 13.0 (3d, pbns, lam) 


Fig. 14: Dynamic Pressure Contour-Case 5 (Pa) 


Fig. 15: Dynamic Pressure Contour-Case 6 (Pa) 


VI. Total Pressure Contours 



g 2 96e + 00 
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1 32e+00 
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J 


Contours of Total Pressure (pascal) 

May 15.2016 


ANSYS FLUENT 13.0 (3d. pbns. lam) 


Fig. 16: Total Pressure Contour-Case 1 (Pa) 



Contours of Total Pressure (pascal) 

May 16. 2016 


ANSYS FLUENT 13.0 (3d. pbns, lam) 


Fig. 17: Total Pressure Contour-Case 2 (Pa) 
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Contours of Total Pressure (pascal) 

May 16. 2016 


ANSYS FLUENT 13.0 (3d, pbns, lam) 


Fig. 18: Total Pressure Contour-Case 3 (Pa) 



Contours of Total Pressure (pascal) Jun 25. 2016 

ANSYS FLUENT 13.0 (3d. pbns, lam) 


Fig. 20: Total Pressure Contour-Case 5 (Pa) 
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Contours of Total Pressure (pascal) 

May 16.2016 


ANSYS FLUENT 13.0 (3d, pbns. lam) 


Fig. 19: Total Pressure Contour-Case 4 (Pa) 



Fig. 21: Total Pressure Contour-Case 6 (Pa) 


VII. Calculation of Head Loss Coefficients 

The head losses in the individual branches can be calculated using the following formula [1]: 

, (Pt 1,2,3- Pt Inlet) 
k = — j 

(1) 

Where; 

P T i , 2, 3 Total Pressure in branches 1, 2 and 3 
Pt inlet Total Pressure in Inlet Pipe 
Vt inlet Reference flow velocity at Inlet 
p -> Density of air 

The velocity magnitudes and pressure values obtained from the fluent analysis are used to carry out the calculations of head 
losses for all the branches of each trifurcation case. The reference inlet flow velocity (Vi n i e t = 3 m/s) [2] and density of air (p 
= 1.225 kg/m 3 ) are constant for all the calculations. 

The head loss coefficients have been calculated and their values have been tabulated in the table XII. 
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The above formula yields negative values of head loss coefficients for all the branches of trifurcations. However, the non- 
dimensional coefficients (k) can be called energy change coefficients rather than head loss coefficients whenever branching 
of flows occurs [3]. Thus, the negative sign can be ignored here and head loss coefficients can be considered as energy loss 
coefficients. 


Table 12 

Head Loss Coefficients in the Branch Pipes of Asymmetric Trifurctions 


Trifurcation Case No. 

Head (Energy) Loss Coefficients (k) 

Branch- 1 

Branch-2 

Branch-3 

1 

2.60 

2.20 

2.65 

2 

2.63 

2.23 

2.72 

3 

2.69 

2.26 

2.73 

4 

2.67 

2.25 

2.53 

5 

2.70 

2.26 

2.64 

6 

2.72 

2.30 

2.69 


VIII. Discussion of Results 

The velocity magnitude contours, dynamic pressure contours and the total pressure contours for all the six cases are as shown 
in the figures 4 to 9, 10 to 15 and 16 to 21, respectively. It has been observed that the velocity magnitudes, dynamic and total 
pressures in the branch pipes decrease with the increase in the trifurcation angles for all the six cases. And also maximum 
values have been observed in the branch-2. 

From the analysis of the contours for all the above six cases, it can be seen that the values of velocity magnitude, dynamic 
and total pressures are maximum in the central branch (branch-2) of the trifurcation. The distribution of velocities, pressures 
and separation of flow in the other two branches of the trifurcation mainly depends upon the turbulence at pipe trifurcation 
junction, angle of trifurcation, and diameter ratio [4]. 

The values of energy change coefficients have been calculated for all the six cases of trifurcations and are tabulated as shown 
in the table XII. It can be observed that the branch-2 head loss coefficients are smaller compared to that of the other two 
branches. This is because there is only change in the pipe area and more energy dissipation that is taking place is because of 
the viscous friction at the wall, while the side branches suffer a directional flow change (secondary flow) along with the 
cross-sectional variation [1]. 


IX. Conclusion 

Energy (head) loss analysis has been carried out and volumetric flow rates, velocity magnitudes, and dynamic & total 
pressures have been determined for all the six cases of asymmetric trifurcations. It is seen that the fluid flow rate, velocity 
magnitude and dynamic & total pressures are more in the branch-2 compared to the other two branches. Smaller values of 
head losses have been obtained for branch-2 in all the cases. This is because there is only change in the pipe area and more 
energy dissipation that is taking place is because of the viscous friction at the wall, while the side branches suffer a 
directional flow change (secondary flow) along with the cross-sectional variation [1]. The turbulence at pipe trifurcation 
junction, angle of trifurcation, and diameter ratio are mainly responsible for the losses and separation of flow [4]. 
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